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ABSTRACT 


This thesis investigates the application of adaptive 
filtering at the NUWES test ranges. Two adaptive algorithms, 
least-mean-squares and recursive-least-squares are studied. 
To facilitate the investigation, a model of the test ranges 
was developed. This model accounts for spherical spreading 
and linear attenuation of the propagated acoustic signals as 
well as the effects of doppler shift, multipath, and finite 
propagation delay time. After describing the model, the 
adaptive filtering algorithms are briefly explained. Then, 
two schemes of adaptive filtering, adaptive noise cancellation 
and adaptive line enhancement, are applied to the model. 
Simulation results of the noise cancellation and line 


enhancement schemes are presented for several scenarios. 
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I. INTRODUCTION 


A. THE NUWES TEST RANGES 

The NUWES ranges are testing and evaluation facilities for 
the Navy. They are comprised of a series of computer linked 
underwater hydrophone arrays. The hydrophone arrays, in turn, 
provide tracking of test vehicles launched from known 
locations. The tracking signals emitted by the test vehicles 
are received by the hydrophones and processed by the 
computers. The recovered signal provides identification and 
telemetry information. 

As part of the scenario, a broadband countermeasure 
produces noise that disrupts the tracking signals. This 
interference reduces the test ranges' ability to recover 
vehicle identification and telemetry information. 

In an effort to lessen the effect of the countermeasure on 
the telemetry, some form of Signal recovery is necessary. 
This thesis investigates one type of signal recovery called 


adaptive filtering. 


B. PROCEDURE 

In this thesis, an appropriate model of the NUWES test 
ranges is first developed. This model accounts for signal 
propagation loss, wideband noise sources, propagation signal 


delay, signal multipath, doppler shift, and the frequency 


response of the hydrophones. Next, a brief explanation of the 
recursive algorithms is given after which, the algorithms are 
developed and applied to the model. Finally, detailed 
Simulation results of the tests conducted on several NUWES 


test range scenarios are presented. 


C. ORGANIZATION 

The thesis contains five chapters and three appendices. 
Chapter II develops the physical model of the NUWES test range 
as well as the signal model used in the simulation. Chapter 
III gives a brief explanation of the LMS and RLS algorithms. 
Chapter IV develops the application of adaptive algorithms to 
the model developed in Chapter MII. Conclusions’ and 
recommendations on the use of the algorithms at the NUWES test 
ranges are presented in Chapter V. Three appendices are also 
provided. Appendix A contains additional results that 
supplement the discussion in Chapter IV. Appendix B furnishes 
a flow chart of the major programs used in the simulation. 
Appendix C contains a listing of the programs mentioned in 


Appendix B. 


II. SIMULATION MODEL FOR THE NUWES TEST RANGE 


A. BASIC SETUP 

The simulation model is based upon the NUWES test ranges. 
The ranges consist of numerous hydrophone arrays arranged on 
the ocean floor. A typical layout of one of the ranges is 
shown in Figure 1. The circular areas indicate the regions 
within which the hydrophone arrays can reliably receive a 
Signal. The area covered by the reception circles surrounding 
the hydrophones delineates the test range. 

An acoustic source emitting a 75 kHz BPSK tracking signal 
is attached to each test vehicle. The tracking signal is 
received by the hydrophones and processed by station computers 
to produce tracking data. 

In a noiseless environment, all test vehicles can easily 
be tracked; however, the addition of countermeasure noise 
interferes with the tracking signal and corrupts the tracking 
davea.. 

To recover the corrupted tracking data, some type of 
filtering is needed. This thesis investigates the use of 
adaptive filtering to enhance the tracking of test vehicles in 
the presence of countermeasure noise. 

To aid the investigation and analysis of signal recovery, 


a model of the NUWES test ranges was developed. This model 





Figure 1. A hydrophone arran gement typical of the NUWES test 
ranges. 


& 


accounts for: 


¢ Tracking signals 

¢- Wideband countermeasure noise 

- Spherical spreading 

- Linear attenuation 

¢ Multipath effects from surface reflection 


* Doppler shift for fast moving targets 


B. HYDROPHONE GRIDS 

The model contains four subarrays arranged in a pattern 
typical of the NUWES test ranges. With these four subarrays 
and their attached hydrophones, any scenario in which 
countermeasure noise is a factor can be examined. 

Figure 2 displays a mesh plot of the subarray layout. 
Each hydrophone subarray has a maximum effective range of 1500 
meters. The range of coverage of each array is depicted by 
the raised circular portion of the mesh plot. The distance 
between hydrophone subarrays is 2500 meters. The hydrophone 
subarrays themselves are depicted by the small protrusion in 
the center of each circle. 

The hydrophones and associated preamplifiers at each array 
are physical devices with finite frequency response 
characteristics. In the simulation, the frequency response 
of the hydrophone/preamp combination was modeled by a twelfth 


order Butterworth bandpass filter. The filter's magnitude and 
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Figure 2. Mesh plot of the subarray layout. 


phase response are illustrated in Figure 3. These 


characteristics closely match those of the actual system. 


C. SIGNAL MODEL 

The tracking signal is modeled as a 75 kHz BPSK signal 
with a code length of 47 bits. Each bit lasts 93 us, 
resulting in a signal bandwidth of approximately 20 kHz. The 
countermeasure is a broadband jammer modeled as white noise. 
The signals are sampled at 300 kHz at the output of each 


hydrophone. The tracking signal is 


x) =ACOsS (2m £36) (2-1) 


Where £.1s the carrier frequency (75 kHz) and A is a square 
wave of fixed magnitude representing the binary code, whose 
Slgn is positive for a 'one' bit, and negative for a ‘zero' 
Bue. 
1. Signal Attenuation 

Signals lose power as they propagate through the 
water. The transmission loss may be considered to be the sum 
of two types of losses: spreading and attenuation. Spreading 
loss accounts for the weakening of a signal as it spreads out 
from a source. Attenuation loss accounts for absorption and 
scattering. Each of these losses in Signal power is 
proportional to the distance traveled before reception at a 


hydrophone. (Urick, 1983, pp. 99-103) 
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Magnitude and phase response of the hydrophone 


Figure 3. 
model. 


Spherical spreading is a type of spreading loss in 
which the signal power is reduced as it propagates unbounded 
from a source in an unbounded medium. The Signal is assumed 
to propagate radially, and the intensity of the signal 
decreases as the square of the range. As the power radiates 
outward it must cover increasingly larger areas; this reduces 
the power at any one point. The transmission loss due to 


spreading is 


Loss = 10 logr?® = 20 logr (2-2) 
where r is the range from source to receiver. (Urick, 1983, 
lo. 200-101) 
Unlike spreading loss, absorption loss’' varies 


logarithmically with range. At frequencies above 10 kHz, the 
attenuation of underwater signals due to sound absorption 
becomes significant. Figure 4 (Clay and Medwin, 1977) shows 
the attenuation profile of underwater acoustic signals in dB/m 
versus frequency in kHz. Thus, the absorption loss for a 75 


kHz Signal is 


a=0.04r aB (2-3) 


where r is the range. 
Figure 5 depicts the power attenuation versus distance 


for a 75 kHz acoustic signal propagating through the test 
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Figure 5. 


range. Both types of losses, spherical spreading and 
attenuation, are accounted for in the illustration. 
2. Propagation Delay 

In salt water, sound travels at about 1500 meters per 
second (Urick, 1983, p.322). This relatively slow speed can 
cause significant delay times between transmission and 
reception of a signal. If a tracking signal source is at the 
outer edge of a hydrophone coverage area, the delay would be 
ds Per. Propagation delays of this order can destroy any 
correlation that might exist between signals received at two 
different hydrophones. Thus, since correlation between 
Signals received at two different locations is a critical 
feature of adaptive filtering schemes, propagation delay 
becomes an important consideration. The delay time is given 


by 


Te (2-4) 





where T is in seconds and r is the range traveled in meters. 
3. Multipath 

Signals arriving at a receiver via multiple paths can 

cause constructive or destructive interference. Since the 

hydrophones are mounted on the ocean bottom, only the multiple 

paths caused by surface reflection of the signals is 


considered. These multipath effects are taken into account 


eZ 


when generating both the tracking and countermeasure noise 
Signals. A delayed and attenuated version of the noise signal 
is added to the original noise signal. However, the tracking 
Signal is transmitted at discrete time intervals and the 
reflected signal is typically received in an interval that 
does not overlap with the time interval in which the direct 
Signal is received. Thus, after the reception of the direct 
path signal, the subseguently received reflected tracking 
Signal can be identified and safely ignored. 
4. Doppler Effect 

The freguency of a received signal depends on the 
relative motion between source and receiver. The relationship 
between the transmitted frequency and the received frequency 


werpgiven by (Urick, 1983, p. 322) 





f! = i) (2-5) 


where c 1s the sound speed and v is the relative speed of the 


vehicle. The frequency shift is 


Af = +0.69Hz/ (knot) -(kHz) . (2-6) 


With a maximum speed of 50 knots (about 26 m/s) for the 
vehicle, the associated doppler shift at 75 kHz is + 2.6 kHz. 
The sign of the shift depends on whether the signal source is 


moving toward or away from the hydrophone. 
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III. ADAPTIVE FILTERING 

Adaptive filtering provides a method of recovering a 
desired signal in additive noise when the signal statistics 
are either unknown or are slowly varying with time. If the 
first and second order signal statistics were known, then a 
fixed optimal filter such as the Wiener filter could be used. 
However, without complete knowledge of the signal statistics, 
some form of estimation of those statistics is required. Two 
algorithms that help estimate filter coefficients under these 
conditions are least mean squares (LMS) and recursive least 
Squares (RLS). (HaykKin, 1984, p. 2) 

The Wiener-Hopf algorithm provides the basis of these two 
methods. The Wiener-Hopf equation for estimating a finite 
impulse response filter weight vector a is given by (Haykin, 


1984, p. 32) 


a=RyL,, (3gf9 


where R, is the autocorrelation matrix of the input signal and 
Xr, is the cross-correlation vector between the desired signal 
d(n) and the received input signal x(n). Adaptive filtering 
schemes fundamentally try to approximate the weight vector of 


the optimal filter by adjusting it "on the fly". 
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A. LEAST MEAN SQUARES (LMS) ALGORITHM 

The LMS algorithm is based upon the method of steepest 
descent and the knowledge that the mean-square error forms a 
paraboloid in the filter coefficient space. The error e(n) is 
defined as the difference between the desired signal d(n) and 


the filter output y(n) (see Figure 6) 


e(n) = d(n) -y(n) (32) 


where 


y(n) = xT(n)a (323) 


and x(n) is the vector of observations (the received signal). 


Thus, the mean-squared error is 


J = Eje*(n)]. (3-4) 


A plot of J versus a, called the error surface, is a 


paraboloid in shape. Expanding equation (3-4) gives 


Ue Ona leider a. (3-5) 
Differentiating the mean-squared error with respect to the 
weight vector produces the gradient vector V 


+2R 2. (3-6) 


d te. 





The basic adaptive filter model. 


Figure 6. 
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The gradient vector points in the direction of steepest ascent 
on the error surface. Based on the gradient search or 
steepest descent method the filter weight vector is 


recursively updated as follows 


a(n+1) = a(n) +=p[-V(n)] (3-7) 


where - V(n) points the filter weight correction toward the 
Minimum point on the error surface and wp is the step size 
parameter. Successive corrections in the direction of the 
negative gradient vector eventually minimize the mean-squared 
error and lead to an optimum filter weight vector. (Haykin, 


moe DD. 93=—102) 


B. RECURSIVE LEAST SQUARES (RLS) ALGORITHM 

While the LMS algorithm is based on ae stochastic 
formulation, the RLS algorithm is a finite data formulation 
based method. It 1s computationally more expensive than the 
IMS algorithm, but it is known to converge to the optimal 
solution significantly faster than the LMS algorithm. (Haykin, 
1984, p. 149) 

The sum of the squared errors between the desired signal 


d(n) and the filter output y(n) is given by 
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y= Dei = Vi(a-a's,). (3-8) 


a= ea 


Differentiating with respect to the weight vector and setting 


the resulting equation equal to zero gives 


OG Anns: 7 if 
ae 22 es ae X= 0. (3-9) 


Solving equation (3-9) leads to 


Re eae or (3-10) 


where 


R, =o (%, 82) (3-10a) 
<2 


Ww 


is the input autocorrelation matrix and 


ec: | (3-10b) 


is the cross-correlation vector between the filter input and 


the desired signal. Now, Rand r can be written as 
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ne4 


1 = Ob (ESHPS) ee Ss eraser (3-11) 
rs 
and 
nea) 
xr, = >) («,4;) +44, = a Bese ar Co) 


igi Uine Pee Fase (3-13a) 
SubSeleucing “Rivas. = ro. and then adding and SUbtrace rng 


x.x'.a , on the right-hand side of (3-13a) yields 


Kna, = RQ, *%,(In~ X44): oa 


n 


The term in parenthesis on the right-hand side of equation (3- 
13b) is a close approximation to the defined error and is 


denoted by e Simplifying (3=-13b) and solving for the 


(aint) * 


weight vector gives 


(3-14) 


Equation (3-14) provides a recursive formula with which to 


update the weight vector. However, evaluating R,' is 


og 


computationally expensive. The matrix inversion lemna, 
commonly used to recursively update the inverse correlation 


matrix, 1S given by 


Ree (7?) Raa, 


Re =7Raee = : ~ (3-15) 
1 exeaa) Ree, x (72) 
For notational simplicity, define 
Ri, X(N : 
‘, 6. Sa ee ie (3-16) 
dt eae era Ta) 
Substituting (2-16) into (3-14) and) (3-15), yielae 
a, 7 Be * Kn 2 cyin-1) (350 
and 
R, = Races ke Coa (3=ie> 


respectively. Equations (3-17) and (3-18) provide recursive 
formulas for updating the weight vector. (Haykin, 1986, pp. 


385-387) 


C. RLS VERSUS LMS 
There are some basic differences between the RLS and LMS 
algorithms. Both of these algorithms recursively adjust the 


welght vector to minimize the mean-squared error. However, 
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the algorithms differ in the way they update their respective 
weight vectors. 

The LMS algorithm updates its weight vector based on the 
product of the error signal, the filter input, and the step- 
size parameter. Furthermore, the step-size parameter, along 
with the input signal, determines the speed of convergence for 
the weight vector. The choice of the parameter is constrained 


such that 





Or 1 x (3-19) 


wneues, is the largest eigenvalue in the correlation matrix 


x 


R,- The product of the scalar step-size parameter and the 
gradient vector produces a vector with which the filter weight 
vector is adjusted. (Haykin, 1984, pp. 100-103) 

While the LMS filter weights are adjusted using a scalar 
transformation of the observations, the RLS filter weights are 
adjusted using a matrix transformation. Also, the LMS 
adjustment vector contains information available from the 
current iteration only, while the matrix used in adjusting the 
RLS filter weights is an estimate of the autocorrelation 
matrix for all of the past data. Therefore, the RLS filter 
weights undergo a higher order adjustment then the LMS filter 


weights, causing the RLS method to converge faster. (Haykin, 


1984, pp. 148-149). 
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The superior performance of the RLS algorithm, however, is 
attained at the expense of increased computational complexity. 
The number of multiplications per iteration required for the 
RLS algorithm is on the order of 3M(3 + M)/2 where M is the 
filter order. In contrast, the multiplications per iteration 
required for the LMS algorithm is on the order of 2M + 1. 


(Hayking 1984, pee 149) 
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IV. ADAPTIVE FILTERS FOR THE NUWES TEST RANGES 

This chapter discusses two adaptive filtering schemes, the 
adaptive noise canceler and the adaptive line enhancer. Each 
scheme is developed from a simple lossless model to a complete 
model for the NUWES test range. Simulation results are 
presented for various scenarios in this development. Note 
that the term "noise" used throughout the chapter refers to 
the countermeasure signal, and the signal-to-noise ratio (SNR) 


is defined as 








2 
SNR = 10LOG,, estexiert) (4-1) 
9) 
where E[(x’?(n)] 1S Signal power, and o? is the noise 


(countermeasure) power. 


A. ADAPTIVE NOISE CANCELER 

The adaptive noise canceler consists of two inputs. The 
first input, the primary input, contains the desired signal 
corrupted by noise. The second input, or reference input, is 
derived from a source at which the tracking signal is weak, 
but contains noise which is highly correlated with the noise 
in the primary input. The reference Signal is adaptively 


filtered to maximize its correlation with the primary signal, 
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then subtracted from the primary input to cancel the 
countermeasure noise. (Widrow, 1975, 1692-1693) 
1. Noise Cancellation of an Ideal BPSK Signal 

The adaptive noise canceler configuration for a 75 kHz 
BPSK signal in broadband white noise has the form depicted in 
Figure 7. Since the signal and noise have comparable 
propagation delay times, zero propagation delay is assumed for 
both the signal and countermeasure noise in the initial 
Simulation. The signal received at the tracking array is the 
tracking signal corrupted by countermeasure noise. At the 
same time, the reference array receives only the 
countermeasure noise. Once received, the signals are passed 
to the adaptive noise canceler as the primary and reference 
Signals, respectively. Figure 8 displays the time and 
frequency plots of the received and recovered signals. As 
seen in the signal spectrum plots, the noise that was present 
in the input signal is attenuated in the recovered signal. 

2. Doppler Shift with ANC 

The ANC scheme shown in Figure 7 can be used to 
process the doppler shifted BPSK signals. The noise canceler 
does not remove the doppler frequency shift, but removes the 
noise from the frequency shifted tracking signal. Figure 9 


illustrates the recovery of a doppler shifted BPSK signal. 
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3. Multipath Effects on ANC 

Signals transmitted from an underwater source can 
propagate to the receiver via multiple paths. In this thesis, 
it is assumed that the surface reflected signal and the direct 
path signal comprise the multipath signal. The reflected 
Signal is a delayed, sign inverted, and attenuated version of 
the direct path signal. Both the countermeasure noise signal 
and the tracking signal travel via these multiple paths. The 
reflected countermeasure noise increases the overall noise 
variance, the amount of which depends on the path lengths and 
the associated attenuation. 

Unlike the countermeasure noise signal, the tracking 
Signal is sent in discrete bursts with quiet time in between 
bursts. As such, the reflected tracking signal, when 
received, does not overlap in time with the direct path signal 
and can be disregarded by the NUWES tracking station. Figure 
10 depicts a noise cancellation scheme incorporating the 
multipath situation. Since the noise is white, the multipath 
noise in the desired signal is uncorrelated with noise in the 
reference. The received signal spectrum in Figure 11 
illustrates the addition of correlated and uncorrelated noise 
to the desired signal. In the recovered signal spectrum, the 


noise uncorrelated with the reference is still present. 
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10 GB). 


4. Signal Propagation Delay 

The ANC functions on the premise that some correlation 
exists between the noise in the primary Signal and the noise 
in the reference signal. However, signal propagation delay 
destroys the correlation between broadband noise signals. 
Therefore, to allow the noise canceler to function with delay, 
the noise correlation must be restored. To reestablish the 
noise correlation, either the primary or the reference signal 
must be appropriately delayed. For realtime tracking, the 
delay must be applied to the reference signal. Delaying the 
primary signal would not permit realtime processing. 
Therefore, if delaying the reference noise, the countermeasure 
noise propagation time to the reference receiver must be 
shorter than the propagation time to the primary receiver. 
Thus, some form of local reference is required for the 
countermeasure noise. The local reference can be obtained by 
placing a receiver near the countermeasure noise source. 

If realtime tracking were not a consideration then one 
of the nearby non-tracking hydrophone subarrays with a large 
noise signal could be used as the reference receiver. In this 
case, the first signal received, primary or reference, would 
be delayed to achieve maximum correlation. 

5. The Overall Model for ANC 
Figure 12 illustrates the overall simulation scheme 


incorporating delays, multipath, and correlation. The 
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tracking signal and the countermeasure noise propagate with 
delay to the tracking array. In addition, a reflected, 
attenuated, delayed version of the countermeasure is received 
at the tracking array. The sum of these three signals 
comprise the corrupted tracking signal, which is in turn 
passed through a bandpass filter representing the tracking 
array. This filtered signal represents the primary input to 
the adaptive filter. The reference signal varies from the 
primary signal because it contains no tracking signal. 

The correlator depicted in Figure 12 adjusts the delay 
of the signal received at the reference array to maximize the 
correlation between the primary input signal and the reference 
input. The error between the filter output and the primary 
input signal becomes the recovered tracking signal. 

Figure 13 shows the output of the system at a signal- 
to-noise ratio (SNR) of -10 QB. (Additional graphs at 
different values of SNR are included in Appendix A.) This 
scenario presumes that correlation between the noise in the 
primary signal and the reference signal exists only between 
the direct path noise. Therefore, the recovered tracking 


Signal contains the uncorrelated, surface-reflected noise. 


B. ADAPTIVE LINE ENHANCER 
The Adaptive Line Enhancer (ALE), illustrated in Figure 
14, differs from the adaptive noise canceler in that the 


reference signal x(n) is derived from the input signal d(n) 
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uSing a delay. 

The input signal is a narrowband signal corrupted by white 
noise. Since white noise has a narrow autocorrelation 
function (theoretically an impulse), the delay used in 
producing the reference signal removes the correlation between 
the noise components of the primary and reference signals. 
However, the narrowband signal has a wider autocorrelation 
function and retains its correlation despite the delay. 

The output of the adaptive filter is the correlated 
portion of the primary and reference signals. This output is 
subtracted from the primary input signal to minimize the error 
Signal. 

The ALE scheme used to recover the tracking signal is 
shown in Figure 15. The input to the ALE is the sum of the 
BPSK tracking signal and the broadband countermeasure noise. 
However, the filtering effect of the hydrophone reduces the 
bandwidth of the countermeasure noise to about 20 kHz. 
Therefore the noise in not broadband, nor is the tracking 
narrowband. In fact, the noise bandwidth is equal to the 
bandwidth of the tracking signal. Since the ALE is effective 
only with a narrowband signal in broadband noise, the ALE will 
not be able to recover the tracking signal. Figure 16 shows 
the output of the ALE when trying to recover the tracking 
Signal. 

A whitening filter inserted after the hydrophone will not 


improve the recovery of the tracking Signal. Since after the 
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Figure 16. The received and recovered signals for the basic 
adaptive line enhancer scheme (SNR = -10 GB). 
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hydrophone, both signals have the same bandwidth, the 
whitening filter broadens both signals, and the ALE will not 


be able to recover the tracking signal. 


C. PROBABILITY OF ERROR COMPARISONS 

Monte Carlo simulations were run using the complete signal 
models developed above. The simulations used signal-to-noise 
Batios of 0, <-5, -10, -15, and -20 QB. The recovered 
Signals were heterodyned with a coherent source then low pass 
filtered to remove the higher frequency components. Once the 
bits were recovered, they were compared with the original bit 
pattern to determine the number of bits in error. 

At each SNR level, for both ANC and ALE schemes, the 
process was repeated 100 times. The tracking signals were 
represented by a random bit pattern of 47 bits. Thus, at each 
SNR, there were 4700 bits. The total number of bits in error 
throughout the 100 trials were accumulated. Figures 17 and 18 
display the percentage of bits correctly recovered as a 
function of SNR. On both figures, the 'x's mark individual 
recovery percentages for each of the 100 trials while the 
solid line is the overall percentage of bits correctly 
recovered. Also, for comparison, probability of error results 
were obtained for the case in which no adaptive filter was 
used. From this case, the ALE and the ANC schemes are 


compared as shown in Figure 19. 
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Figure 17. Probability of error for the adaptive 
cancellation scheme. 
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enhancer scheme. 


41 








N 
L} 
vy 
\ 
7 cf 
S 3 
= — 
5) = 
= © oc 
S|; a 
ra \ 
jes 
E SZ 
= \ Te 
© \ 
: \ 
5 \ s 
ne = 
SS 
. 
= | 
} | 
'@) | St 
O 
Fa \o 
_. i" 
‘0 
! ! S 
—N 
= a) nN WV of wr CS Von 
a S oC S = S = 
=) = = = 
PI9A0d9y 1U9919g 
Figure 19. Probability of error comparison for the 


filtering schemes versus the ‘no filter! case. 


42 


The adaptive noise canceler consistently provided about 4 
dB of Signal gain above the 'no filter' case. At 0O dB SNR the 
probability of an error is about zero percent. As the SNR is 
decreased to -20 dB, the probability of error increased to 
about 23 percent. 

On the other hand, the adaptive line enhancer attenuated 
rather than recovered the tracking signal. The performance of 
the ALE for this type of signal was worse than if no filter 
had been used. This poor performance was due to the broadband 


nature of the tracking signal. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY OF RESULTS 

This thesis investigated adaptive filtering for the 
recovery of vehicle tracking signals in the presence of 
countermeasure noise. A model of the NUWES test ranges was 
developed to which the adaptive filtering schemes were 
applied. This model included a 75 kHz BPSK signal to simulate 
the tracking signal and white noise to represent the 
countermeasure. These signals were subjected to propagation 


losses, propagation delays, multipath effects, and doppler 


shifts. 
Two adaptive filtering methods were examined: adaptive 
noise cancellation and adaptive line enhancing. The first 


scheme, adaptive noise cancellation, employs two inputs. The 
first input, the primary input, contains the tracking signal 
corrupted with countermeasure noise. The second input, the 
reference input, contains countermeaSure noise only. Since 
the noise in the reference input is correlated with the noise 
in the primary input, the ANC cancels the noise, thereby 
improving the received tracking signal. 

In the simulation, a reference receiver was placed near 
the countermeasure noise source to ensure that the propagation 


time for the reference input was shorter than that for the 


44 


primary input. Therefore, the delay used in aligning the 
Signals was applied to the reference and not to the primary 
input, allowing realtime processing of the tracking signal. 

Results of a Monte Carlo simulation are depicted in Figure 
19. The horizontal distance between the ANC line and the '‘'no 
filter' line is a measure of the processing gain. The 
adaptive noise canceler provided about 4 dB of processing gain 
across the 0 to -20 dB SNR range. The results show that the 
correlator used in the noise canceler was unable to correlate 
the surface reflection signals. 

The second filtering scheme investigated, the adaptive 
line enhancer, recovers narrowband signals from corrupting 
broadband noise. However, the BPSK tracking signal used by 
the NUWES test stations has a bandwidth of over 20 kHz. 
Furthermore, the filtering effect of the hydrophone diminishes 
the bandwidth of the countermeasure noise to 20 kHz. The 
Similarity of these two bandwidths prevents the ALE from 
recovering the tracking signal. Figure 19 shows that the line 
enhancer fails to provide any signal gain. In fact, signal 


recovery worsened when using the line enhancer. 


B. RECOMMENDATIONS 

For the model developed in this thesis, the adaptive noise 
canceler provides significant noise attenuation. However, to 
implement the filter, minor modifications to existing range 


equipment must be made. To track in realtime, some form of 
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local countermeasure noise receiver must be used to furnish 
the reference input. If realtime tracking were not a 
consideration, then one of the nearby non-tracking hydrophone 
Subarrays with a large noise signal could be used as the 
reference receiver. In this case, delay is applied to 
whichever signal is received first, primary or reference, to 
achieve maximum correlation. 

The adaptive line enhancer would be easier to implement at 
the NUWES test ranges as it requires only one receiving 
hydrophone array. However, an adaptive line enhancer is 
effective only with a narrowband signal in wideband noise 
(Haykin, 1984, pp. 18-19, Widrow, 1975, pp. 1711-1712). Thus, 
recovery of the signal in its present form is not feasible 
with the line enhancer. However, if the signal could be 
spectrally compressed into a narrowband form, adaptive line 


enhancement could be a viable filtering scheme. 
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ADDITIONAL RESULTS FOR THE ANC 


APPENDIX A. 
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APPENDIX B. PROGRAM FLOW 


A. ENVIRONMENT PROGRAMS 


CETUE: Sets the parameters describing the environment 
| ----- SS Creates the environment 
Ea BOTM: Places the hydrophones 
| ----- SCENE: Calculates SNR's and delays 
: ----- RZ: Creates the signal sources 
| aeeae TRAKARAY: Determines the 


| tracking array 


| meee DISTANCE: Calculates distances 
| to the arrays 


: -<<—<-- LVLS: Computes SNR's at the 
| arrays 


> 


B. FILTERING PROGRAMS 


DOIT: Sets the parameters describing the filtering 


| ----- FILT: Simulates the filter 
| 


Pee DESIRED: Produces the signals used in the 
! filter's primary input 


| 
| 
| 
| 
| 
| 
| 
| 
| | eee RCVD 
| 
| 
| 
| 
| 
| 
| 


Ee REFERENCE: Produces the reference input 


| 
| 
| 
| 
| | 
| a FILTER 
| 


| ----- ADAPT: Filters using either RLS or LMS 


54 


APPENDIX C. PROGRAM LISTINGS 


A. PROGRAMS THAT CREATE THE SCENARIO 


1. Setup 


%* This is the main program that sets the scenarios of the 


% NUWES tests 
clear 

% variables: 
cmdpth=200; 
sdpth=100; 
emcol=35; 
cmrow=20; 


scol=25; 
srow=35; 


° 
© 


ZHREKKKKKKEKKKKKKKES 


the depth of the countermeasure 
the depth of the signal 


oP oO 


the starting position for the 
countermeasure 


oO oO 


oo 


the starting position for the signal 


call the functions *kRKRKKKKKKKKKKK 


[z,zs,zcem,sdelay,cmdelay,SNRS,TA,RA]J=nuwes... 
(EMmaptn, sapth-cmeol- cmrow;,scol) Ssrow); 


oe 


Z 

ZS 

zcom 
sdelay 


cmdelay 


O° dP of cf oO O© © oO oO dO oO oP a oP AO cl oP oP OP O° 


1 - snrdat 
2- snrrt 
3 - snrdr 
4 =- snrrr 


Returned variables: 


= the grid bottom with no signals for 
mesh plotting 

= the signal source vrbl for mesh 
plotting 

7 the interference vrbl for mesh 

plotting 

matrix of the delays for the signal 

source 

= matrix of the delays for the 
countermeasure 


xxzk*k Note ***x* the delays returned are relative. The 
shortest delay in each matrix has been subtracted from all 
elements in the corresponding array. This reduces the 
actual size of the data set required to model the system 
without changing the model. 


SNRS is a matrix of the four following values in order 


SNR for the tracking array direct path 
SNR for the tracking array reflected path 
SNR for the reference array direct path 
SNR for the reference array reflected path 


SD 


HP dO HO 


JO WO BO dO oO oO of O1O CH 1° OO 


xxxk The above order is followed in all vectors used *** 
TA = the chosen tracking array number 
RA = the chosen reference array number 


col ==> 0 oeeeweereee ee 8) 
QOW sooo nee ---- == 

0 | Al | 

| | 

A3 
| | 
A2 

| | 

—_ Ad | 
TQ), SS SS Se 


2. Nuwes 


function [(z,zs,zcm,sdelay,cmdelay,SNRS,TA,RA}=nuwes... 


0 oO Bd oP oO GO CO co BO ce dO ol oO GO oP oP oP of GO Of OP OO CLO oo GO ao CH oO 


(cmdpth,sdpth,cmcol,cmrow,scol,srow) 


variables: 
cmdpth = the depth of the countermeasure 
sdpth = the depth of the signal 


the starting position for the 
countermeasure 


Emecol*> eEmrov 


scol; srow = the starting position for the signal 

Z = the grid bottom with no signals for 
mesh plotting 

ZS = the signal source vrbl for mesh 
plotting 

zZcom = the interference vrbl for mesh 
DLoOcLtLing 

sdelay = matrix of the delays for the signal 
source 

cmdelay = matrix of the delays for the 
countermeasure 


*kkk Note *** the delays returned are relative. The 
shortest delay in each matrix has been subtracted from all 
elements in the corresponding array. This reduces the 
actual size of the data set required to model the system 
without changing the model. 


SNRS is a matrix of the four following values in order 

1 - snrdt SNR for the tracking array direct path 

2 - snrrt SNR for the tracking array reflected path 
3 - snrdr SNR for the reference array direct path 

4 - snrrr SNR for the reference array reflected path 
*kxkk The above order is followed in all vectors used *** 
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% TA = 

% RA = 

% 

% Sole — 0 
% row --- 
% 0 ! 

: | 

% 

2 | 

: 

: | 

% : 

% 10, == 
araydpth=400; ss 
soundspd=1500; % 

M=(-30 30]; % 
sdb=186; % 
cmdb=sdb+20; % 
Xis=('square') ; 


load bottom 


% botm creates: 

% X = the 
% Y = the 
% s1,2,3,4 = the 
% 

% xmax = the 
% dir 
% ymax = the 
% dir 
% dx = the 
% the 
% dy = the 
% the 
% Z = the 


the chosen tracking array number 
the chosen reference array number 


ARRAY LAYOUT 


A3 


A2 


the depth of the array 

speed of sound in water is 1500 m/s 
the mesh orientation 

the signal level in dB re uPa 

the countermeasure level in dB re uPa 


[X,Y,S1,S2,83,S4,xmax, ymax,dx,dy,z]=botm; 


xX portion of the mesh grid 
y portion of the mesh grid 
position of the arrays in meters as 


measured from the top left corner (0,0) 


max number of mesh intervals in the x 
man number of mesh intervals in the y 


spacing in meters between intervals in 

x direction 

spacing in meters between intervals in 

y direction 

botm grid with radius circles and arrays 


{SNRS,cmdelay,sdelay,TA,RA,zs,zcm]=scene... 


(cmdpth,sdpth, 


araydpth,soundspd,sdb,cmdb,cmcol,... 


Cmrow ;SCol{srow,X%, Y ;dx,aGy;,sl,s2,s3,s4); 


3. Botm 
function [{X,Y,s1,s2,s3,s4,xmax, ymax,dax,dy,z]=botm 
% X = the x portion of the mesh grid 
% Y = the y portion of the mesh grid 
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s1,2,3,4 = the position of the arrays in meters as 
measured from the top left corner (0,0) 

the max number of mesh intervals in the x dir 

the man number of mesh intervals in the y dir 

the spacing in meters between intervals in the 

x direction 

the spacing in meters between intervals in the 

y direction 

the botm grid with radius circles and arrays 


xmax 
ymax 
ax 


dy 


Of BW oO © oO oO oO oO AO 


Z 


this subprogram will return the bottom configuration of 
the NUWES test range. The model of the bottom will 

contain only four sensors. This is sufficient to model 
all desired cases of signal vs countermeasure position. 


create the bottom grid of four offset sensor arrays with 
a 1500 m detection radius. Each sensor is 2500 meters 
apart measured center to center. 

The upper left hand corner is (0,0). 


Of oO lO Oo oO GO GO oO oO CO 


The positions of the sensors in (x,y) coords with x being 
the horizontal axis are: ( in meters) 


of o? 


$1=(1500 1500];s2=[(4000 1500];S3=[(2750 3665] 7s4=[(5250 3665]; 


$3 


Of © oP oP co olf oP oP oP 


% 
% Now, the mesh for the bottom will be created 
% The grid is sampled every 10 meters to decrease the 
% reguired points 
xmax = 70; % the maximum number of samples in the x 
% direction 
ymax = 70; % the maximum number of samples in the y 
* direction 


range=7000; # the range length 
adx=range/xmax; % the distance between x sample points 
dy=range/ymax; % the distance between y sample points 
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¢ Set up the x-y grid 


xg=1:xmax;X=ones(xg) '*xg; 


yg=1: ymax; Y=yg'*ones(yg) ; 
Z=zeros (xg) ; 


$ 
$ now put in the rings for the sensors 


for x=1:xmax 
for y=1: ymax 
1f ((x*dx-sl1(1))*2 + (y*dy=-sl1(2))*2) <= 1500%2 
Z (xy) 5 


elseif ((x*dx-s2(1))%2 + (y*dy-s2(2))*2) <= 1500%2 


zZ(X,y)=5; 


elseif ((x*dx-s3(1))*2 + (y*dy-s3(2))*2) <= 1500%2 


Z(X,Y)=5; 


elseif ((x*dx-s4(1))%*2 + (y*dy-s4(2))*2) <= 1500%2 


2(xX,Y)=5; 
else 
Z2(x,y)=0; 
end 
end 
end 


$ put in the array markers 
z2(s1(1)/dx, s1(2)/dy)=10; 
z(Ss2(1)/dax, s2(2)/dy)=10; 
Zess(1)/ax, s3(2)/dy)=10; 
Zesa(1)/dx, s4(2)/ady)=10; 


4. Scene 
function [SNRS,cmdelay,sdelay,TA,RA,Z2s,zcm]=scene... 


(cmdpth,sdpth,araydpth,soundspd,sdb,cmdb,.. 
cmcol,cmrow,scol,srow,X,Y,dx,dy,s1,s2,83,84); 


% variables: 

$ SNRS is a matrix of the four following values in order: 
$ 1 - snrdt = SNR for the tracking array direct path 
$ 2- snrrt = SNR for the tracking array reflected path 
$ 3 - snrdr = SNR for the reference array direct path 
% 4 - snrrr = SNR for the reference array reflected 

% path 

$ cmdelay = matrix of the delays for the 

$ countermeasure 

$ sdelay = matrix of the delays for the signal 

% source 

% TA = the chosen tracking array number 

% RA = the chosen reference array number 


52 


ARRAY LAYOUT 


A3 


j | 
! | 
! | 
= Gee | 
! | 
| { 


AO AP AO AO dO AO AP CO AO AP ONO 


SESEEEEETEEETEEEETEEETEEETETEEES 
SCECECESCEEEEEEEE ESTEE ETE ETELEEEES 


EEE 
43%$%% 
and 


cm, zcm]=RZ(xX,Y¥7ax, ay, cmcol ,cmrow ema); 

create a noise source and place it within the bottom 
mesh. Then find ranges to noise source and the dB re uPa 
level at each range usSing spherical spreading and linear 
attenuation. 


of Wo oO HOP dO oO HM CO 


CECEEEEEEEEEELEEEEELEETESE SESS 
for signal source %%ESE¥EEEESF 
SEEEEEEEEEEEETEETEEFETEETEE TSS 
ay aed ane eee 

create a signal source and place it within the bottom 

* mesh. Then find ranges to noise source and the GB re uPa 
+ level at each range using spherical spreading and linear 
% attenuation. 


O° o\ oO 
of oP oP 


>< ole 


Crs A olf oO 
J of ol of? 
W 
~~ 
N 


CECCCECECCCECCECCEEEECESCCESESEEEEEEEEEEEEEEEEEFLESEEEELEEESESESS 
%$%3%S%S%S%SSSS¥ESEES% Find the tracking array tEVESEEESEVESESSSSS 
EECECEECECECECECEECESEEEEEESEETSEEEEEEEEEEGEEEEEEEEEESSESTESS 
%* The tracking array is the array closest to the signal 
source 

% 

(TA,RA,dsdt,dcmdt,dsdr,dcmdrjJ=... 


trakaray(sl,s2,s3,s4,Rs,Rcm,Zs,zcom,dax,dy) ; 
% find the array closest to the signal source and the 
 non-tracking array with the lowest SJR 


% A = matrix of array positions 

$ TA = the chosen tracking array number 

$ RA aa the chosen reference array number 

% dasdt = direct distance to tracking array from 
% the signal 

% dcmdt = direct distance to tracking array from 
% the countermeasure 

% dsdr = direct distance to reference array from 
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the signal 
direct dist to reference array from 
the countermeasure 


of off AP ol? 
QO, 
Q 
J 
QO, 
Ky 
ll 


EELELELEEEEEELEEESEEEEEEEEEEEE EEE EEE EEE EEE EEEEEEEEEEELEEEES 
%$%%4%%¥E¥EE¥EE¥EESE Find all of the distances and times tEEEEEEEES 
ESSSEELEESSSELSESEEESSSEEELEEEEELE EEE ELE EEE EE EEE EEL EE EE EEE ESS 


% general form 

$ (dtd,dtr,drd,drrj=distance(d2t,d2r,araydpth,dpth) 
$ dtd = direct distance to tracking array 

$ Ger = reflected distance to tracking array 
% ard = direct distance to reference array 

% anc = reflected distance to reference array 
% 

% Gaze = flat distance to tracking array 

% a2r = flat distance to reference array 

% araydpth = depth of the array 

% dpch = depth of the signal source 

% NOTE: All depths and distances are in meters 


oe 


find distances and times for the countermeasure 


femeta ;cmatr,cmdrdad, cmarr |=... 

distance (dcemdt,dcmdr,araydpth,cmdpth) ; 
cmdist=[cmdtd,cmdtr,cmdrd,cmdrr]; 
cmdelay=cmdist/soundspd; 


[sdtd,sdtr,sdrd,sdrr]=distance (dsdt,dsdr,araydpth,sdpth) ; 
sdist=(sdtd,sdtr,sdrd,sdrr]; 
sdelay=sdist/soundspd; 


% decrease the delays such that the shortest delay is zero 
cmdelay=cmdelay - min(cmdelay) ; 
sdelay=sdelay - min(sdelay) ; 


SSESELSELTEESESESELELEEEE EEE EEE EEE EEE EEE EEE EE ELE EE EEE ELEES 
SELEEEEEEEEEES compute the SNR levels FESEEEEEEEELELELEES 
ESEEEEEEEEESE ES at the arrays EEELESEEEELEEEEEEES 
EESEEEEEEELEEELEEEL EEE EE EEE EEE EEE LEE ELE EEL EEEE EEE EEE ELEES 
% 
general form: 

function [snrd,snrr)=lvls(sdb,cmdb,dsd,dcmd,dsr,dcmr) 


lvls computes the signal to noise levels for the 
reflected and the direct paths to the array for which the 
%$ distances are given. 


% 
$ 
% 
$ 
$ 
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% variables: 

% snrd - SNR for the direct path to the array 

% snrr = SNR for the reflected path to the array 
% 

% sdb - the dB level of the signal source 

% cmdb = the dB level of the countermeasure source 
% asd = distance from the signal direct 

% dcmd = distance from the countermeasure direct 
% dsr = distance from the signal reflected 

% dcmr = distance from the countermeasure 

% reflected 

% 


* find the SNR for the tracking array 
(snrdt,snrrtj=lvls (sdb, cmdb,sdtd,cmdtd,sdtr,cmdtr) ; 


% find the SNR's for the reference array 
(snrdr,snrrrj=lvls(sdb,cmdb,sdrd,cmdrd,sdrr,cmdrr) ; 


SNRS=[snrdt,snrrt,snrdr,snrrr]; 


B. PROGRAMS THAT PRODUCE THE SIGNALS AND DO THE FILTERING 


1. Doit.m 


“4 elle ie 

% variables 

% SNR Signal to Noise Ratio 

$ £) filter level 

% D the delay for the adaptive line enhancer 

% old old = 1 then use previously calculated vrlbs 
% old = 0 create new vrbls 

% seed the seed value for the noise generation 

% pass pass = 1 then pass just the filter coeff and 
% inverse corr matrix to next trial (saved in 
% filtcoeff. ie save filtcoeff a Rinv) 

% code 47 ones and zeros as desired 

$ filttype filter type 

% O = lms 

% 1 = rls 

% blocks the number of blocks of data produced 

$ D the delay for the adaptive line enhancer 

% D= 0 no ALE 

$ D> O use the ALE 


$clear 

%+SNR=; $ set by start.m 
$fl=; $ set by start.m 
$D=; $ set by start.m 
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oo 


$seed=; set by start.m 


old=0; Stk Olde ——l chem Vebis- mat must be Ehe variables 
% from the previously run case. 

pass=0; 

$code=(ones(1,47)]; 

eoae—-{] 1°) 0 1°0 02 2 7°02 0°07 1.101001 1 =... 
Or OO 1 tery O 2-00, brit 0 1.0: 0-1 a yO 1); 

filttype=1; 

blocks=1; 


(sn,dn,yn,en,J,seed,code,xn,f,fs,a,h,t,bits,Rinv,f1,SNR,.. 
tbase, bitsbase,blocks]=filt(SNR,f1,D,old,seed,pass,code,... 
pa letype, blocks) ; 


Save vrbls % vrbls.mat is used when old=1 
% sn the signal from the torpedo 
% t the time index 
% if the waveform frequency 
% fs the sampling freq 
% J the jammer signal 
% J2 the delayed jammer Signal (set delay below) 
% sd the propagated torpedo signal 
% Ss} the propagated jammer signal 
% dn the desired received signal sd+sj 
% bc,ac the butterworth filter coeffs 
% yn the filter output 
% en the desired signal - the filter output 
% xn the reference input 
% a the adaptive filter coefficients 
% h the computed wiener filter coefficients 
% bits the bit index for the x-axis 
% Rinv the correlation matrix for the adaptive 
% filter 
2. Filt 
f u n c 1a 1 fe) n 


(sn,dn,yn,en,J,seed,code,xn,f,fs,a,h,t,bits,Rinv,f1,SNR,... 
tbhase,bitsbase,blocks]=filt(SNR,f1,D,old,seed,pass,code,... 
filttype,blocks) 

this program runs the filter routine 


calls: desired produces the desired signal 
and assorted variables needed 
throughout the program 


ref produces the reference signal 


OP oP oP oO AP AO AO oO 
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adapt does the adaptive filtering and 
produces the output 


taps calculates the wiener filter tap 
weights for the noise assuming 
the noise is stationary for the % 
sample period 


0 AO A WO AO OO 


EESEEEEEEEETEEEEEET ESTEE TEEETEL ESTELLE ET EEEETEE TEE EE LEE EEEESS 
$4%%%%%% PRINT OR NO PRINT AS THE PROGRAM PROGRESSES %%%%%% 


if pl = 1 then plot graphs. match pl in the 


pl=0; $ 
: % next loop down 


ESEECECEEEEESESSECECETEEEEEEEEEEEEELEEEEEEEEEEEEEEEEESETETESS 


$S%ESSESEEEEEESEEEESE CALLING DESIRED %E%EtEVEVEVELELESEEESS 
% DESIRED will produce the waveforms used throughout the 

% program 
ESEECEEEEESECESEECEESTESCEEEEEEELEVEEEEEEEEEELELEVEVESESTETETS 
if pl== 

msg='entering desired.m' 

end 


$$%SS%%SS%%S6666% Load the old vrbls here tEttttEEFtEEESEEF 


% Old variables are loaded when passing values from one 

% program run to the next program. Use old variables when 
% you want to expand the data set w/o increasing the sample 
% size. 


O then create new values 
1 then load the last set of vrbls 


if old==1 % old 
% old 


load vrbls; 
pl=1;0ld=1; 


tmax=t(length(t)); % get the max time in the last 
% iteration 
t=tbase + tmax; % correct time for number of the 
% last iteration 
bitmax=bits(length(bits) ); %* get the max bit number 
bits=bitsbase + bitmax; % correct bit number for 
% iteration 
[rw,co)=size(a); 
a(1,:)=a(rw-fl,:); % this sets the filter 
% coeffs to the final 
% values from the last use. 
else % create new values 


(sn,t,f,fs,J3,32,sd,sj,dn,bits,bc,ac,tbase,bitsbase]=... 
desired(SNR,code,blocks, seed) ; 
end 
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$%4%%%S%%E%S% end of 1if-then FEBEEEELEEES 
* desired creates: 
% sn the signal from the torpedo 
% t the time index 
% fe the waveform frequency 
% fs the sampling freg 
% J the jammer signal 
% J2 the delayed jammer signal (set delay in 
% desired.m) 
% sd the propagated torpedo signal 
% sj the propagated jammer signal 
$ dn the desired received signal sd+sj 
% after the bandpass filter 
% (1f you want no BPF then go into desired.m 
% and comment 
% out the filter line and activate the line 
% below it) 
% Also, you can add a reflected jammer into dn 
% bits the bit numbers for the x axis 
% bc,ac butterworth filter coeffs 
% tbase the time base 
% bitsbase the base scale for the bits 
if pl==1 
msg='leaving desired.m' 
end 


%$%%$%% LEAVING DESIRED %%%%%%%%%%%%%%% 
% CEEETEEELEESEEESE ESE ESEEEFESECSS% 


oo 
ov 
ove 
ove 
ove 
ove 
ovo 


$%%S%SSSESSES6ESESSSSS6S6S6%S6%6% CALLING REFERENCE %%%%%% 


Tf pl== 
msg='entering ref.m' 
end 


if old == 
1f D== 
[xnJ=ref(J,bc,ac) ; % change J to J2 when time delayed 
% change J2 to J when no time delay 
% set the delay in desired.m 
% change J to J+J2 when adding a 
% 


reflected jammer 


HreaakKKKKKKKKKKKKKKKKKKKKKKKK KKK KKK KKK KKK KKK KKK 


% The following linea are for the adaptive line enhancer 
LekkkkkkkkkkkkkKeeekKekKeekkkeekakkeekkkaekkekkekekkkeekkKeekkKeKek 


elseif D>0 
% {dn)=hdyne(dn,t,f) ; % mult the received signal 
%* by the local osc and put 
%* thru LPF 
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delay=D; % the amount of delay for 
the ALE 


xn=[zeros(1,delay) dn(1:length(dn)-delay) }; 


end 
PEST CCCCCCTCCCOCCC OCS S SS SS CC CSCS SSS TCT CCS SSCS SSCS SL Se ST 


% end of the adaptive line enhancer variable adjustment 
FRAAKKKKKAKKKKKKKREKKKKKKEKKKEEKEKEKEEKEEKEKKEKKKKKKKKKKKKKK 


end 


% reference creates: 


% xn the reference signal for the adaptive filter 
if pl== 

msg='leaving ref.m' 

end 

ESESESCEESEEEESEESEEEEEEEEEEE% LEAVING REFERENCE %2%%3%4%%%SESS% 
EECECECESECECECEECECEEECESEEEFEEEEEEFESEEEFESEFESESESTEESTS% 
EEEESEEEEEESEESSESEEEEEESEEE CALLING ADAPT FESESEEEETESEEEFEFSS 
if pl==1 

msg='calling adapt.m' 

end 


if old== 

* the program will enter this loop only if 

* new data is being used and the filter is adapting from 
tap weights of all zeros (initial conditions). If the 
* program does not enter this loop, then the filter will 
adapt to the data using and improving upon the old tap 
% weights. 


start the filter coeffs with 
all zeros 

setup of the inverse 
correlation matx 


a(1,:)=zeros(1,f1); 


Rinv=200*eye(fl,fl); 


Of AO WO OW 


if pass == pass = 1 means that the filter 
coefficients are passed to the 
next array but the noise and 
other signals might have changed 


load filtcoeff retrieves ‘a' and 'Rinv'! 


dO do AO W AO ot 


[rw,co]=size(a); make the initial weights for the 
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%* filter equal to the 
a(l,:)=a(rw-fl1l,:); % final values from the last use. 


end 
end 
EEEEESSESSESSE EEE TC SEFESEEE EEL EE TES EEL EEEE EEE SELES 
$EEE%SS%ESEE Choose an adaption method FEETEETELESES 


if filttype == 


(en, yn,aj=adaptlms(dn,xn,fl,a); lms method 
elseif filttype == 

fen,yn,a,Rinv]=adaptrls(dn,xn,fl,a,Rinv) ; rls method 
end 


% Creates: 


$ en the error between the desired output 

% and the filter output. This is our 

% filtered signal 

% yn the output of the adaptive filter 

% a the filter coefficients. Each row 

% is the filter coefficients as time 

% progresses. 

% Rinv the inverse of the correlation matx. This 
% variable is passed to avoid the large 

% variations at the start of a new iteration. 


save filtcoeff a Rinv % variables used to pass parameters 
%¢ between arrays or between trials 


end 
E%%SEEEEEESSSS¥SES¥SSESSX¥SES¥E LEAVING ADAPT EEEEEEEESESSESSSESSS 
ER CESSES ESCES SESE EES LESSEE ESE SEE EEEE EEE EE EEE EEE SEE EEE ES SEES ESS 


SCESETESSEEEELESEELEESEEEE ENTERING TAPS EEEEEEEEEEEELES ESS 
$ Computes the wiener filter tap weights 


if pl== 
msg='entering taps' 
end 
*(hj=taps(sn,xn,dn,t,fl); %h = tap weights for wiener 
% filter assuming that the 
$ noise is stationary 
% in the sampling period 
if pl== 
msg='leaving taps' 
end 


SEEEEEEEEEEEEEETEEEEEEEEE LEAVING TAPS EEEEEEEELESEETEF 
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3. Desired 


function f{sn,t,f,fs,0,J2,Sd,S),an- bites be, ac, tbase, 
bitsbase ]=desired(SNR,code,blocks, seed) 

% assembles the desired signal 

to produce the signal from the torpedo 

to produce the noise from the jammer 

to produce the received signal at the 


% calls: wave 

% jam 

% rcva 

% 

% bandpa 
% 

% filter 
% 

% Creates: sn 
% NE 
% pe 
% fs 
% af 
% J2 
% sd 
% sj 
% dn 
% bc 


array 


ss to create the filter coeffs modeling the 
hydrophone 
bandpass filter to the signals 


applies 


the 
the 
the 
the 
the 
the 


Signal from the torpedo 

time index 

waveform frequency 

sampling freq 

jammer signal 

delayed jammer signal (set delay 
below) 
the propagated torpedo signal 
the propagated jammer signal 
the desired received signal sd+sj 

Jac the butterworth filter coeffs 


{[sn,t,fs,f,bits]=wavecode(code,blocks) ; 
produces the signal from the torpedo 


provides 

% sn 
% s 
% fs 
% bits 


%SNR in dab used 
tbase=t; 


bitsbase=bits; 


K=(10* (-SNR/20)) 


the signal 

the time index 
the sampling freq 

the bit numbers for the x axis 


to scale the noise 


+ the base set of time used for adding to 
% iterations 
%* the base set of bits used for adding to 
% iterations 


VSQGEe (2). 


J=K*ejam((t t],seed)/2; 


J2=J3 (1317:1317+1 


J=3(1:1316); 


315); 


[bc,ac}]=bandpass(f,fs); 


dP dP dP dP dP dP AP dP 
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produces the noise from the 
jammer 

produces the noise for the 
reference with a delay 


produces the bandpass 
filter coeffs 

bec = the numerator coeffs 
ac = the denominator coeffs 


path=1; 

[sdj=revd(path,sn); % produces received signal at the 
array 

path=2; 

[sj J}=rcevd(path,J); % produces received jammer at the 

array 

path=3; 

[s}2}=rcevd(path,J2); * produces received reflected jammer 

* at the array 

SESEEEEELEELELEEEELEETELELE LEFT SESS 

lgnal (comment out three of the four) 

% (remember to change the reference signal in filt) 

adn=filter(bc,ac,sd+s)j); % use this line for BPF 

dn=filter(bc,ac,sd+sj+sj2); use this line for BPF with 


reflected jammer 


° 
x) 
S 


Ke) 
ry 
O 
©, oo? 
Cc 


% 
% 
dn=sd+sj; %* use this line for no BPF 
$dn=sd+sj+sj2; %* use this for reflected 

%* Jammer w no BPF 
CCCCECEECEEESCCEEEEESEESEESTETEEFEELEEEEETE ETE ETE TE EL ELEETSETS 


% path = propagation path 
% sa = the propagated torpedo signal 
$s) = propagated jammer signal 
% dn = the combined signal + noise 
% after the bandpass filter 
4. Ref 


function (xnj=ref(J,bc,ac) 
* assembles the reference signal 


$ calls: revd to produce the received signal at the 
% array 

% 

Seer LOrS: none 

% 

* Creates: xn: the signal received at the reference 

% 

path=3; % chooses the propagation path 


{xn)=rcvd(path,J) ; 
xn=filter(bc,ac,xn); 


% produces received signal at the 
reference 

% path = propagation path 

J = original jammer signal 

% xn = the combined signal + 

% noise after the bandpass 

% filter 
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5. Adaptlms 


function [en,yn,a]=adaptlms(dn,xn,fl,a) 


¢ this function will adapt the filter parameters using the 
% Widson-Hopf algorithm (LMS) 

% 

* dn = the desired signal 

¢ xn = the reference signal 

* f1 = the filter order 

% u = acceleration parameter 

% a = the initial filter coefficients 

%* get the max length of the data vectors for loop counter 


nend=length (xn) ; 

if length(dn) > nend 
nend=length (dn) ; 

end 


%¢ pad the data vectors with enough zeros for the loop 
xn=(xn zeros(1,nend-length(xn)+f1)]; 


dn=(dn zeros(1,nend-length(dn)+f1l)]; 


Leekkk LMS METHOD xkxeakake 

¢ enter the loop. the loop adjusts the filter coefficients 
* with each iteration and computes error: dn-yn=en. 
u=.0005;sigsq=0; 


for n=l:nend 


yn(n) = a(n,:) * xn(n:nt+fl-1)!'; % yn for each time step 

en(n) = dn(n) - yn(n); 

delofa= -2 * en(n)*xn(n:nt+fl-1) ; + delta of a used to 
adjust 


% filter coeffs in LMS 

$%%%% The forgetting factor t%t%%%% 

Sigsq=.9*sigsg + xn(n:n+f1l-1)*xn(n:n+fl-1) '; 

un=u/sigsqd; 

a(nt+1,:)= a(n,:) - u*delofa; % The next filter coeffs 
end 

6. Adaptrlis 

function [en,yn,a,Rinv]=adaptrls(dn,xn,fl,a,Rinv) 


%¢ this function will adapt the filter parameters using the 
%* Recursive Least Squares algorithm 
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the desired signal 

the reference signal 

the filter order 

the error signal (our output) 
the adaptive filter tap weights 


OP AO AP AO of OO 
rh 
~ 

tou ue wu 


% get the max length of the data vectors for loop counter 
nend=length(xn) ; 
if length(dn) > nend 
nend=length (dn) ; 
end 
+ pad the data vectors with enough zeros for the loop 


xn=[(xn zeros(1,nend-length(xn)+f1)); 
dn=[(dn zeros(1,nend-length(dn)+f1)); 


3 zkkkkk RLS METHOD xexkxeekae 
% enter the loop. the loop adjusts the filter coefficients 
% with each iteration and computes error: dn-yn=en. 


for n=l1:nend 


ran) = 
en(n) = dn(n) - yn(n);3 


k = ((Rinv * xn(n:n+fl-1)') / (1 + xn(nintfl-1)* ... 
Rinvy*xn (Wins fl=1) 7) 3 
Rinv = Rinv - k*xn(n:n+f1l-1) *Rinv; 


anti, +) =) (a(m,:)'+ kK * (dn(n)=-xni(nen+fl-1)*a(n-_:).") 


end 


71 


ais) “A ea te ae % yn for each time step 
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